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Recently, much attention has been paid on modified gravities, and specifically on their possible role in the 
accelerated expansion of the universe, a fact widely accepted by a great part of the scientific community, and strongly 
supported by the observations. In particular, f(R) theory is the easiest and most popular modification of General 
Relativity (GR), specially because of its simplicity and ghosts absence (for a recent review of f(R) gravity, see 
Ref. [l|). In addition, the reconstruction of f(R) theories that are capable of reproducing the dark energy epoch, 
and even the inflationary phase, is generally straightforward in comparison with other modifications of gravity (see 
Refs. @tJ])- In this sense, f(R) gravity is able of reproducing ACDM epoch (see Ref. Q), or it can directly mimic a 
cosmological constant at the current era, unifying the entire cosmological history[4j. Moreover, some f(R) gravities, 
so-called viable models, can pass the local gravitational tests and reproduce a realistic cosmological evolution (see 
Ref. |3|). As well viable f(R) gravity may pass the matter instability p, a large perturbation that may occur in the 
study of celestial body solution. Nevertheless, viable f(R) models are known to reproduce a kind of future singularity 
(see Ref. p}), a problem that can be circumvented in the scalar-tensor representation of f(R) gravity, |j|. Other 
modified gravities (like Gauss-Bonnet ones) which are also capable of reproducing the dark energy epoch [9( have 
been also studied. 

However, much less attention has been paid to more complex theories, specifically to those theories that contain 
a non-standard coupling of curvature with energy-momentum tensor. Recently, class of modified gravity theories 
in which the gravitational action contains a general function f(R,T), where R and T denote the curvature and 
the trace of the energy- momentum tensor respectively, has been proposed (see Ref. |lQ(). Cosmological evolution 
of such theories has been studied, including the reconstruction of cosmological solutions or the presence of future 
singularities (see Ref. [ll[). Nevertheless, such strong coupling of the curvature and the trace T implies the violation 
of the usual continuity equation, an issue that can be solved by an appropriate function f(R,T), as shown for 
the first time in Ref. [12J- Nonetheless, this function gives rise to cosmological perturbations evolutions for the 
sub-Hubble modes which are in contradiction with the standard GR results, eventually leading to singularities of 
matter perturbations (see Ref. |12| | ) and to appearance of extra force in geodesic equation if conservation law is violated. 

In the present work, an extension of f(R,T) gravity is proposed. The action described by a generic function 
f{R,T,R^ v T 111 ') is considered, where R is the scalar curvature, T is the trace of the energy-momentum tensor, and 
R^ v is the Ricci tensor. Then, the general FRLW field equations are derived in the presence of R^T^ V coupling 
terms. Several cosmological solutions are studied, where it is found that generally the matter sector does not behave 
as in GR, since the divergence of the field equations does not lead to the usual continuity equation. Nevertheless, 
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this issue can be solved by assuming a particular f(R 1 T,R tlv T^ u )^ where due to the extra degree of freedom, the 
usual cosmological evolution for a perfect fluid can be imposed. In addition, a pure de Sitter solution is discussed. It 
can be realized even in the presence of a non-constant fluid. The occurrence of ACDM universe is also analyzed: the 
corresponding gravitational action f(R,T,R^T^) is reconstructed, giving GR with a cosmological constant just as 
particular solution. Finally, the matter instability, a perturbation that may be induced in the interior of celestial 
body solution, is discussed for such class of models. 

Note that physical motivation for the theory under consideration comes from covariant Horava-like gravity with 
dynamical breaking of Lorentz invariance. Actually, such power-counting renormalizable covariant gravity [141 repre- 
sents simplest, power-law version of /(i?,T, i? pi/ T A " y ) theory. One can expect that such theories may provide deeper 
connection between modified gravity and Horava-Lifshitz theory. 

II. f{R,T,R^T^) GRAVITY 

Let us start from the general action for f(R, T, R^ U T^ V ) gravity, 

S = i / ^ x ^~3 f( R > T > R^T^) + J d i x^~g C m = S G + S m . (1) 

Here, k 2 is gravitational coupling constant, R is the Ricci scalar and T represents the trace of the energy-momentum 
tensor, T — Tff, while C m is the matter lagrangian. As usually the energy- momentum tensor is defined as, 

Then, by assuming a matter action that only depends on the metric but does not on its first derivatives, the energy- 
momentum tensor yields, 

T"" = g^C m + 2^ . (3) 

Hence, field equations are obtained by varying the action (TTJ) with respect to the metric g^ v , 

SS= 2^f dAx t^^)-^' T ' R ^ TtlU ) + V=9 (f* SR + fr5T + fpSR^T^)] + J d^xS {^gC m ) . (4) 

Here P = R flv T iiV has been defined for convenience and in order to simplify the expressions that are coming below, 
whereas the subscript means variation with respect to P, R and T. Those terms corresponding to the variations of 
the metric determinant . the Ricci scalar and the trace of the energy-momentum tensor have been already found in 
the literature (see Ref. |10l |) 

V^5 = - 2^9 9^9^ . (5) 

f R SR + f T 5T = [R^f R + (g^D - V^/ v )f R + (T^ + 6 M „) f T ] SgT ■ (6) 

where the second term in the rhs of ([6]) has been integrated by parts, whereas the tensor 0^„ is defined as 
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e - - <f*i£ - ^ - ^ - 2 ^i^> ■ ( 7 ) 



The terms corresponding to the variation of Rf lv T ,lv yield 

f P 5{R^T^) = f P {SR^T"" + R^ST"") . (8) 

The variation of the Ricci tensor is given by 

«V = VaST°„ - V V ST% , (9) 

where 

Sr% = \g aX (Sg MXi u + Sg v x ;fi - Sg^, x ) . (10) 



The variation of the first term in ([5} with respect to the metric yields 
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whereas the variation over the second term in ((5J) is given by, 
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Note that for a regular /(i?, T, R lll/ T til ') function, in the absence of any kind of matter, the f(R) gravity equations are 
recovered. The corresponding properties and solutions well studied in the literature on f(R) gravity arc also satisfied 
by f(R,T,R lljU T^ v ) theories in classical vacuum (for a recent review on f(R) theories, see []J). Moreover, here we 
are interesting to study the behavior of this kind of theories for flat FLRW metrics, which are expressed in comoving 
coordinates by the line element 



ds 2 = -dt 2 + a{tf (dr 2 + r 2 d9 2 + r 2 sin 2 9 dip 2 ) 



(13) 



where a(t) is the scale factor. Then, the main issue arises on the content of the Universe through the energy-momentum 
tensor, and consequently on the matter lagrangian C m through the terms R llv T pu and T = T£ and the trace of the 
energy- momentum tensor T = Tff. Since the flat FLRW cosmology (fT3")) is assumed, the content of the universe 
(pressureless matter, radiation,... ) can be well described by a perfect fluid, whose energy-momentum tensor can be 
defined in several ways, which is assumed here to have the form, 



Tpv = PQnv + {p + p)u^u. 



fji^V • 



(14) 



Here p and p are the energy and pressure densities respectively, and u M is the four- velocity of the fluid, which satisfies 



u l _ l u' J 



T, and in comoving coordinates is given by u M = (1 0).By assuming C m — p, the second variation of the 



matter Lagrangian in ([7]) and (J12I) becomes null. Hence, the tensor O^,,, defined in (J7J, yields, 

&fj,u = -2X^1, +p g^ . 
The tensor S M „ defined in (fl"2"|) can be expressed as, 



(15) 



Hpu — —(j/ivP ~ ~zRTfj, v + 2it ^T^) 



(16) 



Hence, the complete set of the field equations is given by, 



R^fR-\g^f+(g^O - V M V„) f R +(T^ + 9^) f T +\ pT^fp + g^V a V fi T^fp)-V a V^ ) fp+E^fp = k 2 T^ , 

(17) 
where recall that P = R lu/ T 1 "' ' . Then, for a particular equation of state (EoS) p = wp, the corresponding FLRW 
equations can be obtained, and its cosmology studied. 

III. COSMOLOGICAL SOLUTIONS 

In this section, some particular cosmological solutions are studied, and the corresponding action is reconstructed. A 
pressureless fluid, where w m = 0, is assumed in this section. From the one side such choice simplifies the calculations, 
and from the other side, a presureless fluid represents a suitable description of the baryonic and cold dark matter 
content of the universe, whereas the dark energy sector is comprised by the modifications of GR. Then, the energy- 
momentum tensor reads T{f — Diag (— p, 0, 0, 0), whereas the trace T and the scalar P give rise to 



T = T^ 



P = RauT^ = RooT 



00 



-3p(H 2 + H) 



(18) 



A. Model f(R, R^T^) = aR + f(R^T^) 

Let us start by considering the action f{R,R iiv T IJ ' v ) = R + f{R^ v T tlu ), where the dependence on the trace T is 
omitted. By assuming the flat FLRW metric (IT5|) . the FLRW equations are 
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At this point, a particular cosmological solution, described by a given Hubble parameter H(t) can be considered. 
The corresponding action /(R^T^) could be reconstructed. In addition, one can note that the usual continuity 
equation is not longer valid. Hence, the time dependence of the energy-density p(t) has to be also obtained from 
the above equations (|19p , which involve more complicated treatment of reconstruction the gravitational theory than 
within the framework of f(R) gravity. 

First of all, we consider cosmological solutions of the type, 

777 

a(t) = a t m ->■ H = — , (20) 

that shall be refereed to as power-law behavior. Within GR, this type of solutions accomplishes the scale factor 
evolution for perfect fluids with a constant EoS, such as dust (m = 2/3) or radiation (m = 1/2) dominated Universe. 
Then, since the FLRW equations (fT9| shall be composed by powers of the time variable, a natural consideration would 
be to assume p(t) and f(P(t)) to be powers of the cosmic time as well, 

p(t) = pot* , / = f3P n = P (R^T^T = [-3pom(m - 1)]" t" ( ^ 2) , (21) 

Then, the FLRW equations JT9) yield 



m 2 2 ^ a [n 2 (<7-2) + 3n(m + l)+m-l] [-3m(m - lp )]" , . 

oa—z — k pot + p — — 1 = u , 

t A 2(m — 1) 

6a(m - l)(3m - 2)m 2 [3(m - l)m + n(2 + m)(l + 3m) + n 2 (cr - 2) (4m - 3) + n 3 (a - 2) 2 ] 

6(m— l)mt 2 6(777—1)777 

x (-3/j rr7(m - l)) n i"( ff - 2 ) = (22) 

One can analyze several cases, depending on the kind of the action considered through the coupling constants choice. 

• In the most general case, where all the coupling constants are non-zero, the only possible solution of the equations 
yields, 

cr = —2 , ?7 = 1/2 . (23) 

Whereas, {a, j3} are expressed in terms of m, 

9m -7 2 „ 2(3m-2)(m-l) , ^_ 

a = ? = -K 2 po , B = i ' \ ' K Jog ■ (24) 

6m(3m 2 -3m + l) '"' ' ( 3 m 2 - 3m + l)^-im{m - 1) V ^ V 7 

In order to ensure a real and physical gravitational action, 777 < 1 according to the expression of B in (J24I) , which 
allows the radiation or matter dominated evolution of GR-type, but forbids any accelerating solution (777 > 1). 

• By considering a = 0, the gravitational action turns out to be / = BP n = B{R^, v T ,J ' v ) n , and from (|22[) . 

n = -?— . (25) 

a — I 

Then, the equations (|2"2"j) give rise to a system of transcendental equations that can be solved numerically for 
{777, a} for each value of 8, leading to a set of solutions. 



Finally, by assuming that the coupling among the gravitational and matter sectors is described solely through 

, basically by setting k 2 — 0. Then, the exponent of the action ([21]) is constrained to, 

2 (26) 
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As above, the equations (|22|) become a system of algebraic equations for the constants {m, a} in terms of {a, (3, }, 
which describe an infinite set of solutions, but confirms the suitability of the action with f(P) = P" = (Rui/T flv ) n 
to reproduce power-law expansion. 

Nevertheless, the reconstruction of the exact expression for the gravitational action is not possible in general, and 
numerical calculations are required, as shown below. In addition, note that the evolution of the energy-density behaves 
differently to the one given in GR. This may not be important for the illustration of the reconstruction procedure, as 
for instance, the solutions considered above, but it may seriously affect the observational constraints when a realistic 
evolution is studied. In this sense, the FLRW equations ([T9l) are two differential equations for the functions p(z) and 
f(R liV T fil ') — f(z), expressed in terms of the redshift. They do not ensure that the solution for the energy-density 
behaves in the proper way as observations suggest, namely p m ex (1 + z) 3 for dust matter, apart from the trivial 
solution f{R llv T^ v ) = — 2A that recovers GR with a cosmological constant, and the usual continuity equation is 
satisfied again. Hence, more general actions should be considered in order to provide a realistic evolution. 

B. Model f(R, T, R^T*") = R + f(R„„T^) + g(T) 

By assuming an action described by f{R,T,R IXtV T^ v ) = R + f{R llv T> iU ) + g(T), the FLRW equations are slightly 
modified in comparison with the above case, 



3H 2 + \ 



-3H 2 - 2H - - f 
2 



f+g- 3Hd t (pf P ) - 3 (zH 2 - ij) p/ P ] - p( K 2 - g T ) = , 

f + g- d tt {pfp) - AHd t (pfp) - (zH 2 + h) pf P ] = . (27) 



Within this section, the above equations are expressed in terms of the redshift 1 + z = —rpc, where 4? = — (1 + z)H-4-, 
whereas P — P(z) and T = T(z). In addition, a new variable is defined x( z ) = pfp which simplifies the equations 
(|27|) by extending the system to a set of three differential equations 

3H 2 + i {/(*) + g(z) - 3 [3H 2 + (1 + z)HH'] X (z) + 3H 2 (l + z) X '(z)} - (V + 9 -^j p{z) = , 
-3H 2 +(l+z)HH'-± {f(z) + g(z) - [3H 2 - (1 + z)HH'] X (z) + [3(1 + z)H 2 - (1 + z) 2 HH'] X '(z) + (1 + zf X "(z)} = , 

xW-^f^O. (28) 

Hence, the system is composed by three equations with four unknown functions {f(z),g(z), X (z),p(z)}. In com- 
parison with the above section, the presence of an additional function g(T) introduces a new degree of freedom that 
allows to impose the usual evolution of the energy-density 

p = poa-^+^=p (l + z)^ + ^ , (29) 

Recall that the reconstruction procedure is restricted to pressureless fluids in this section, so (J2U)) reduces to p ex (l+z) 3 
when dust matter w m = is assumed. Hence, the system of equations (|2"8|) is completely solved by the set of functions 
{f(z),g(z) lX (z)}. In addition, note that the equations (|2"8)l can be reduced to a single fourth order equation in f(z) 
by replacing the second and third equations into the first one. In addition, by assuming ([29]) the expressions (p~8|) 
yield 

T = Tji = -'^$n° n (l + zf , P = R^T"" - R Qa T 00 - -3^f&(l + zf{H 2 - (1 + z)HH')) . (30) 



where $7" = /% , and the subscript o refers to the value at z = zq. Moreover, in order to keep the correct units in 
the action, the functions f(P) and g(T) can be defined more conveniently as, 

f{P) = HlF(y\, 9 { T )=HIg(^\, (31) 



where 



p = _^« To = JJ$n° m . (32) 



First of all, let us consider de Sitter (dS) solutions 

H(z) = H . (33) 

The system of equations (|28|) can be easily reduced to a single equation in terms of F(z) that yields, 

2(1 + zfF^(z) + (1 + z) 2 F< 3 (z) - 10(1 + z)F"{z) + 36F'(z) - 1620 m (l + zf = , (34) 

which is straightforward to solve. The solution is obtained in terms of the redshift, 

F(z) = Ci(l + z) a + (1 + z)' 3 {C 2 cos [w ln(l + z)} + C 3 sin [w ln(l + z)]} + C 4 + 30°, (1 + z) 3 , (35) 

where 

a = -1.327, /3 = 3.414, u = 1.38 , (36) 



whereas Cj are integration constants. Moreover, for the dS solution (|33|) and the energy-density (|29|) . the terms 
P/P a = T/T = (1 + z) 3 , and the function F{R^T^ U ) yields, 

P\° /P\ W3 f„ / W , P\ „ . /W, P\l „ „„ n P 



F(P) = Cl ( - j + ^- j |C 2 cos ^ In - j + C 3 sin ^ In - j | + C< + 3<- . (37) 



Whereas G(T) acquires a similar form 

G(T) = C 1 g) Q/ %g)" /3 {c 2 cos(^ln^) + C 3 sin(|ln|)} + C 4 -3^, (38) 

where (7j are constants that depend on the integration constants C, and the parameters (|3"6")l . Then, the complete 
action can be reconstructed for this illustrative case, where the evolution of the energy -density (|29| has been imposed 
in order to satisfy the usual continuity equation. An interesting remark comes from the fact that in the usual 
approach to modified gravity, as for instance f(R) gravity, pure dS solution (|33l) does not admit the presence of a 
non-constant fluid, since the lhs of the field equations is constant, so the energy-matter sector can not depend on time. 
Nevertheless, this is not the case in f(R, T, P) due to the presence of strong coupling between the energy-momentum 
tensor and the curvature through P = R liV T ,iV , as has been pointed above. 

Let us now consider the reconstruction of a theory capable of reproducing the Hubble parameter described by the 
ACDM model, but in the absence of a cosmological constant, a possibility already explored in f(R) gravity (see Q). 
The following Hubble parameter is considered, 

tf 2 = yP 0a - 3 + |. (39) 

Here, po is the energy-density at the current time (with ao — 1), and A is a constant. By expressing Q39p in terms of 
the redshift 1 + z = -, and the cosmological parameters fi° = „^°, 2 and Q° A — ^4t, the Hubble parameter (l39l) is 
rewritten as 

H 2 = H^ n (l + zf + H 2 n° A . (40) 

Note that the model considered previously, f(R, P) = R + f(P) would give rise to a completely different behavior 
of the energy- density in comparison with (I2TJ1) . except for the trivial solution f(P) = — 2A, where GR is recovered. 



Nevertheless, by assuming the model f(R,T,P) = R + f(P) + g(T), the evolution of the energy-density (J29I) can 
be imposed and the system results in a linear system of differential equations over / and g, as was pointed above. 
Nevertheless, the system (|28jl can not be solved exactly in this case (and in general), but numerical calculations are 
required. By considering the energy- density evolution (|2"5|) and the Hubble parameter (|4"0)) , the functions P(z) and 
T(z) take the form, 



qw4n0 
P(z) = R^T^(z) = _V^(1 + z)- 2fi A ) 



T(z) = - 3 J$n° m (l + zf 



(41) 



Then, the system of equations (|2"8")l can be solved numerically for the Hubble parameter (|40p. which leads to a set of 
functions for F and G, defined in (|31|). For an illustrative propose, some of the solutions are shown in fig. 1, which 
illustrates the kind of gravitational action compatible with the ACDM evolution (|40| . As a particular case, by setting 
all the derivatives of F and G to zero at z = 0, but where F(z = 0) = G(z = 0) = 1, the solution yields to an action 
given by, 



f(R, T, R^T^) = R + f{R^T^) + g(T) = R - 6H n A H% =R-2A 



(42) 



where we have used fl A 



3HS 



This is the trivial case that reduces the gravitational action to GR with a cosmological 



constant. However, there is an infinite number of solutions of the equations (|28[1 which keep the dependence both on 
P as T. 





T/T, (xlO 8 ) 
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Figure 1: The right panel corresponds to the function F(P/Po), whereas the left panel shows G(T/To). Each curve represents a set of 
initial conditions of the differential equations 128P . which have been set up at z = 0, assuming values less than unity for the derivatives 

{F( l '(z), G^ l '(z)} in order to be as close as possible to GR. 

Using same technique one can reconstruct other types of evolution in the theory under discussion. 



IV. MATTER INSTABILITY OF CELESTIAL BODY SOLUTIONS 



This section is devoted to study of so-called matter instability due to the modifications of the gravitational action. 
In the presence of gravitational objects, as the Sun or the Earth, tiny modifications of GR may make the system to 
become unstable in the interior solution, a phenomena called the matter instability (see Ref. [6|). This problem has 
been already analyzed in f(R) gravity, where it has been that matter instability can be avoided by an appropriate 
choice of the function that depends on the Ricci scalar (see Ref. [ij). Furthermore, it has been demonstrated that in 
viable f(R) gravities such instability does not occur as well (for review, see [lj). Nevertheless, the problem becomes 
much more complicated in the case of f(R,T,R liv T fiv ) gravity, since the field equations are very difficult to solve, 
even in the linear approximation, as shown below. 



Let us consider for simplicity f{R,T,R liV T> iV ) = h(R) + f 2 (T) + f 3 (R^T^). For a spherical body solution, 
the field equations have to be solved both in the exterior as well in the interior of the body. In the exterior, where 
T = 0, the field equations (JXTj) reduce to the f(R) gravity ones. The curvature Rq becomes constant and is given 
by the solution of an algebraic equation, being Rq = in GR and in some f(R) gravities. Nevertheless, the interior 
solution may differ substantially from GR, introducing corrections in the curvature that may lead to large instabilities. 

For this analysis, the trace equation is considered, which can be easily obtained by multiplying the held equations 
(IT71) by the metric g^ v 



Rf lR - 2/ + 3D/i R + (T + 9)/ 2t + \aTf 3T + V a V T a ^f 3p + 5/ 3p = « 2 T 



(43) 



where P = R jll JT ilv as defined previously. Then, let us assume that T = Tq describes correctly the interior of the 
celestial body, and according to GR the interior curvature gives 



Rq = -K T , 



(44) 



which is basically the trace equation of GR. In f(R) gravity, the issue arises because the extra terms in the action 
may induce a perturbation on the curvature SR that could grow very fast, giving rise to unstable spherical bodies. 
However, for some specific f(R) functions, most of them called viable f(R gravities, the problem can be circumvented, 
@. In f(R, T, R^nT* 1 ") gravity, the problem becomes even worse, since the strong coupling among the curvature and 
the energy-momentum tensor may prevent To being even a solution of the background equation. Nevertheless, by an 
appropriate f 2 (T) function, the energy-momentum tensor is able to behave as in GR, similarly as shown in the previous 
section in the cosmological context, although the possible presence of a large perturbation SR, and additionally 8P, 
or equivalently Tq5R^ u , may appear. Let us assume that T = Tq = — % is a solution of the background equation 
(|4"3")l . and a small perturbation in the curvature is considered 



R = R Q + SR , P = P + 5P 
Then, the equation (J43)) at hrst linear order yields 



(45) 
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Here recall that j\ — f\ (R) and f 3 — f 3 (P) , and the super- indexes denote derivatives with respect to R and P, such 

as fa' — d l ]\jdR l and fif' = a n f 3 /dP % , whereas the tildes refer to covariant derivatives evaluated in the perturbed 
metric. Note that by setting / 3 = 0, the equation of the stability for f(R) gravity is recovered Ref. |J|. If the 
Hilbert-Einstein action is assumed, /i = R, the solution of the equation (14^1) gives SR — 0, where no instabilities are 
produced, as expected in GR. Nevertheless, the equation ([4T>]) is very complicated to be solved, even to get some 
qualitative information, since it involves not just perturbations of the Ricci scalar SR but also of the Ricci tensor 
SR^v through SP. Since both the Ricci tensor and the Ricci scalar contain second derivatives of the metric, the 
equation (J46I) should be rewritten in terms of the perturbed metric Sg^, which gives rise to a more complicated 
equation. 



Nevertheless, for an illustrative propose, let us consider the model of the previous section, where /i = R, and the 
equation (|4"o) yields 



Tof^asp 



a/3\ ,(3) 



tfVVaTo + [T O 9 a0 + 2Tf) /f V Q P + Vr^ a Tf + /^T Q/3 V 



(2), 



1/3 , A2)rpap, 



VpSP 



So/; 



(2) 



-□ 



(/: 



(2) T 
3 J 



v Q v 



(': 



(2W/9 

-in 



2/ 3 (1) P c 



dP: 



n To/; 



(1) 



2V Q V /3 (T^/ 3 (i: 



h - P)i3 



•Cm — ^To 



SR . 



(47) 



In addition, by considering a constant interior solution T =constant, and consequently Pq ^constant, the equation 
(147) is simplified, 



1. 



f( 2 )r 



+ -To/r^D^P + fPTfVcXpSP + (s /i 2) + 2/ 3 (1) P ) £P 



A (1) - Po/ 



(i) 



I* 



6R 



(48) 



Note that a particular solution of (|4"5|) is given by (5P = <5P = 0, where no instability appears. Nevertheless, in general 
the behavior of the perturbation depends mainly on the term in front of SP in (J48I) , where by imposing positivity on 



such term, Sq/ 



n„ f ( 2 ) . 



2/ 3 



(1) 



Po > 0, the general solution of the equation (|4"5j) would approximately give rise to a damped 



oscillator, avoiding the appearance of a large instability in the interior of an spherical body. This example gives a 
very restricted case, but it provides a qualitative illustration on the reconstruction of viable /(P, T, R^T^) gravities. 
However, a full analysis of the matter instability will probably reveal that a large instability is also common in this 
class of theories, which impose strong restrictions on the viability of the theory, but allows to reconstruct realistic 
f{R,T,R^ v T^ v ) gravities. 



CONCLUSIONS 



In the present work, new version of modified gravity which includes strong coupling of gravitational and matter 
fields, R ilv T >iV has been studied. The physical motivation for such theory comes from the covariant Horava-like 
gravity with dynamical breaking of Lorentz invariance [lj]. In fact, such a covariant power-counting renormalizable 
theory [La] represents simplest power-law F{R,T,R^ V T^ V ) gravity. 

Such modified gravity contains extra terms that allow to reconstruct viable cosmological evolution. In this 
sense, several cosmological solutions have been studied, and corresponding gravitational action is reconstructed. 
Nevertheless, the dynamics of the matter sector, a perfect fluid in the cases studied above, has to be fixed 
previously in order to guarantee the same evolution as in GR, or in other words to satisfy the continuity 
equation. Otherwise, a generic f(R,T,R^ v T^ v ) may give a realistic Hubble parameter H(z), but would give 
rise an anomalous behavior for baryonic, dark matter, and any other perfect fluid present in the equations, as 
shown in the first part of section Mil In addition, this kind of theories allows to reproduce pure de Sitter universe 
in the presence of non-constant fluids what is usually a forbidden solution in GR, and f(R) or Gauss-Bonnet gravities. 

Moreover, the ACDM universe can naturally occur in f(R,T,R tll/ T tiu ) gravity, where the usual evolution for 
dust matter, p ex (1 + z) 3 , is set. This is shown in the second part of section Hill for the model described by the 
Hilbcrt-Einstein action plus corrections accounted by the functions f(R lil /T ttv ) and g{T) has been considered. The 
complexity of the equations does not allow to get an analytical and exact expression, but using numerical results the 
behavior of f{R^, v T^ u ) and g(T) has been obtained, as shown in Fig. 1. 

Finally, the last section has been devoted to the study of matter instability within these theories. The equation of 
the perturbations was obtained. It implies perturbations not just on the Ricci scalar but also on the Ricci tensor. This 
fact makes complicated to get the analytical results since the equation should be reduced to a fourth order differential 
equation on the metric tensor by using the definitions of R and R^ v . Nevertheless, by imposing some restrictions, 
qualitative information has been obtained: it was shown that the reconstruction of f^R^T^) requires additional 
constraints in order to avoid matter instability. Then, some viable models, similarly as those in f(R) gravity, capable 
of reproducing the realistic cosmological evolution could be reconstructed. 

It would be interesting to develop further the study of cosmological dynamics of the theory under discussion with 
attempts to simplify its Lagrangian formulation. This maybe probably achieved via the introduction of Lagrangian 
multipliers 14J and subsequent mapping of the theory to power-counting renormalizable covariant gravity with dynam- 
ical breaking of Lorentz invariance. In its turn this may indicate to nice ultraviolet properties of some /(P, T, R^T^) 
models and its relation with underlying quantum gravity. 
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